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Abstract 



Of 
<D 

\ We discuss the sequence of effective theories needed to understand the quali- 

tative, and quantitative, behavior of real-time correlators (A(t)A(0)) in ultra- 
relativistic plasmas. We analyze in detail the case where A is a gauge-invariant 
ON . conserved current. This case is of interest because it includes a correlation 



recently measured in lattice simulations of classical, hot, SU(2)-Higgs gauge 
theory. We find that simple perturbation theory, free kinetic theory, lin- 
earized kinetic theory, and hydrodynamics are all needed to understand the 
■ correlation for different ranges of time. We emphasize how correlations gener- 

. ically have power-law decays at very large times due to non-linear couplings 

Oh, to long-lived hydrodynamic modes. 

D 

^ ■ I. INTRODUCTION 

> : 

^ ■ There is a rich variety of spatial and temporal scales in hot, weakly coupled, ultra- 

relativistic gauge theories, a selection of which are shown in Table |. Depending on the 
problem of interest, a wide variety of methods can be appropriate for describing the physics 
on different scales: naive perturbation theory, hard thermal loops, dimensional reduction, 
kinetic theory, Vlasov equations, Langevin equations, hydrodynamics, and more. For static 
properties of plasmas, the tower of effective theories relevant to different distance scales is 
well understood 0. These effective theories are all Euclidean quantum field theories, in 3 or 
4 dimensions, and may be systematically constructed using standard renormalization group 
techniques. 

For real-time properties, including the analysis of non-equilibrium response, construct- 
ing an appropriate sequence of effective theories which correctly disentangles dynamics on 
different distance (or time) scales is much less straightforward. Some of the difficulties 
are illustrated by the long confusion over one of the basic time scales of hot non-Abelian 
plasmas: the time scale for non-perturbatively large fluctuations in gauge fields (such as 
those responsible for baryon number violation in electroweak theory) M. To elucidate the 
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real-time dynamics of hot gauge theories, one should understand the tower of effective the- 
ories appropriate for different time scales and, at least in simple examples, calculate the 
appropriate matching between such theories. 



typical particle wavelength 
typical particle separation 
Debye length 
plasma frequency 

non-perturbative magnetic length scale 
mean free path: any-angle scattering 
mean free path: large-angle scattering 
temporal scale for topological transitions 



T 

T -i 



(9T)' 1 
(9T)- 1 

(9 2 T\n) 
(9 4 T\n) 
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(non-Abelian theories) 



(non-Abelian theories) 



TABLE I. Orders of magnitude of various distance and time scales in hot, ultrarelativistic, gauge theory 
for weak coupling. Here, and henceforth, "In" stands for ln(g -2 ), except for the case of the mean free path 
for any-angle scattering in an Abelian gauge theory, where it is a genuine logarithmic infrared divergence. 
(Naturally, c = % = 1.) 

The current interest in hot electroweak baryon number violation has motivated substan- 
tial efforts in numerical simulations of hot, classical, lattice gauge theory [[|-|6|]. For certain 
time and distance scales, it has been argued that classical simulations can adequately re- 
produce the behavior of the underlying quantum plasma of interest (or more generally, that 
certain classical lattice results can be translated into corresponding results for the contin- 
uum quantum theory |7|].) Hence, these simulations are a source of "data" for comparison 
with theoretical pictures of hot gauge theories. In particular, there have been recent mea- 
surements by Tang and Smit || of time-dependent correlations (A(t)A(0)) for several gauge 
invariant operators A. An interesting theoretical goal is then to understand how best to 
describe the physics of, and calculate the behavior of, such time-dependent correlators for 
different ranges of time. 

In order to keep the discussion focussed, we shall specialize to the particular class of 
correlations where the operator A is the spatial part j of a conserved gauge-invariant current 
jp (such as the electric charge current in QED, or fermion number current in QCD). As we 
shall explain below, one of the correlations measured by Tang and Smit falls into this class. 
We shall also focus on the spatial average of such currents, 



as did Tang and Smit. However, the principles and techniques underlying our discussion are 
not limited to this chosen class of examples. As our interest is in the dynamics of hot gauge 
theories, we shall assume that the current j M involves particles which feel a gauge force with 
coupling g. We also assume that there is no chemical potential, so that the total charge 
/ d 3 x j° has zero expectation. 

Fig. |I] qualitatively depicts our basic results for the time dependence of the current- 
current correlator, and the physical pictures best suited to understanding the dynamics at 
different time scales. The purpose of this paper will be to explain this graph. We have 
assumed that the gauge symmetry, and the symmetry corresponding to the current, are 
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FIG. 1. A qualitative picture of the real part of the current-current correlation V(j(t)j(0)) in a hot 
plasma, for the case of bosonic charge carriers. The labels on the axes are only meant to denote orders of 
magnitude, except that the "2T 3 " mark is in fact twice the "T 3 " mark in the small-coupling limit. The 
time axis should be thought of as linear before the break and logarithmic after; the vertical axis should 
be thought of as linear. In is short for ln(g -2 ). The transient shown at t ~ 1/T decays exponentially. If 
the current is carried by fermions instead of scalars, the resulting graph is identical except for the initial 
transient, which starts from minus the "J" 3 " plateau value instead of twice that value. 
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FIG. 2. Similar to fig. [I], but for a classical field theory defined on an infinite spatial lattice with lattice 
spacing a. We have specialized to a current j associated with a non-Abelian symmetry. The Abelian 

case is identical except that <? 18 T n a 8 In 5 / 2 and 11 In 5 / 2 on the left and right axes should be replaced by 
<7 20 T 12 a 9 ln 2 and /3£" 12 ln 2 , respectively. The bottom and left axes are labeled in physical units. The top 
and right axis are labeled in lattice units, where the current j has been multiplied by g 2 as compared to the 
continuum normalization and /3l = 4f z (g 2 Ta)^ 1 . Here, In is short for ln/?L ~ ln[(g 2 Ta) -1 ]. The envelope of 
the transient oscillations decays as i~ 5 / 2 . 
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unbroken. We have also assumed that the current is normalized so that the corresponding 
number operator / d 3 xj° counts each particle as ±0(1) [rather than ±0(g), which would be 
conventional for the electric current]. And we assume that the underlying theory is weakly 
coupled: g 2 {T) < 1. 

As will be discussed later, naive perturbation theory in the underlying quantum field 
theory is sufficient to understand the behavior at early times. At later times, it is appro- 
priate to treat the excitations of the theory as (quasi-)particles with definite position and 
momentum and to apply kinetic theory. At very late times, it is most convenient to focus 
on the long-wavelength, collective, hydrodynamic excitations of the system. As will be dis- 
cussed below, the non-linear coupling of the current to these hydrodynamic excitations gives 
rise to a long-time power-law tail in the decay of the correlation function. 

The behavior of the current- current correlation will turn out to be dominated by excita- 
tions of hard thermal modes — that is, modes with momentum p of order T in the continuum 
theory. As a result, some of the interesting physics of hot plasmas (such as the soft, p ~ g 2 T, 
SU(2)-colored gauge fluctuations responsible for electroweak baryon number violation) will 
not be probed in our leading-order investigation, and the corresponding effective theories 
PJ^Q do not appear in fig. [1]. 

Fig. |2| shows the corresponding result for a classical field theory on an infinite spatial 
lattice. (The finite volume case will be discussed in sec. 0.) We shall later explain the causes 
of the qualitative differences (e.g., the short time oscillations and the different power-law 
decay of the long-time tail) from fig. |l|. The differences between figs. [I] and |2] reflect the fact 
that the current-current correlator is sensitive to the hard thermal modes at all times, even 
in the long time limit. Such modes do not have occupation numbers large compared to one, 
and will never be adequately described by a classical theory. The differing origins of the 
long time tails in figs. [I] and ^] will be discussed in detail in section |TV|. 

In this paper, we focus on a qualitative discussion of the correlator. A quantitative 
treatment of some aspects will be deferred to a companion paper M . Throughout this paper, 
orders of magnitude given for time and distance scales will refer to the case of continuum 
theory unless we explicitly refer to classical lattice theory. In sec. II, we examine the short- 
time current correlations using simple perturbation theory, and discuss the reduction to 
kinetic theory at times large compared to 1/T (continuum) or 1/a (classical lattice). A 
discussion of the leading-order perturbative result may also be found in ref. [|Kj, which 
appeared as our work was being completed. In sec. Ill, we discuss the current correlator 
in the context of kinetic theory, and explain the exponential decay at intermediate times, 
as shown in fig. [I]. Section IV reviews the existence of long-time power-law tails in systems 
with hydrodynamic behavior and shows how they arise in the case of the current-current 
correlator. In particular, we discuss the failure of classical lattice simulations to reproduce 
the long-time behavior of the continuum theory. Sec. V explains how the behavior of the 
current correlation changes qualitatively in finite volume systems. Finally, having (we hope) 
made clear to the reader that classical lattice simulations need not generate the same long 
time behavior as high temperature quantum theories, sec. VI briefly reviews the limited 
class of questions for which real-time classical simulations may be useful. The bulk of the 
paper focuses on currents in theories for which the associated charge is carried by a scalar 
field. (These are the theories for which it is interesting to compare the dynamics of the 
continuum quantum theory with the corresponding classical lattice field theory.) However, 
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all the discussion of continuum quantum dynamics also applies to theories where fermion 
fields carry the global charge of interest. Appendix |A] discusses the small differences between 
results for bosonic and fermionic charge carriers. 

In the remainder of this introduction, we review the correlator measured by Tang and 
Smit. Their interest was in electroweak theory, and they simulated SU(2) gauge theory 
coupled to a doublet Higgs scalar <f>. One of the correlations they measured was of the 
spatial average of the gauge-invariant operator 

W a = ltr([$t(D$) - (D$)t$]r a ) , (1.2) 

where D is the SU(2) covariant (spatial) derivative, r a is a Pauli matrix, and $ is a 2 x 2 
matrix made up of <fi and its charge conjugate: 

$ = (zr 2 0*, 0). (1.3) 

They chose this operator because they were interested in studying the thermal properties 
of gauge bosons. Well into the symmetry-broken phase of the theory, W a is an operator 
that, to leading order, simply creates or destroys massive gauge bosons: W a ~ (4>) 2 A a . 
However, W a has another interpretation: it is the current associated with the global SU(2)r 
custodial symmetry of the gauge-Higgs system. Specifically, the theory is invariant under 
the SU(2) L xSU(2) R symmetry 

$ -> L$R\ (1.4) 

where L and R are independent SU(2) rotations. The SU(2) L symmetry is gauged, and the 
SU(2)r symmetry is an additional global symmetry for which W a is the spatial part of the 
current. 

The SU(2)r symmetry is broken in the low temperature phase of the theory (though the 
vectorial combination of SU(2)r and SU(2)l survives to guarantee equality of the W and 
Z boson masses). In this paper, we shall focus on the hot, symmetric phase of the theory. 
In this case, the operator W a does not simply create and destroy gauge bosons, and it is 
more usefully thought of as a current carried by the scalars. Indeed, W a cannot couple to a 
state composed only of gauge bosons because W a carries an SU(2)r charge [a is an SU(2)r 
adjoint index] while the gauge bosons do not. 

II. PERTURBATION THEORY 

In this and the following sections, we consider Abelian or non-Abelian gauge theories with 
scalar and/or fermionic matter with a conserved, gauge-invariant current j^. The simplest 
example to keep in mind is QED, coupled to either scalars or fermions, where is the electric 
current. For non-Abelian gauge theories, the gauge currents are not gauge-invariant, and 
we do not wish to consider unphysical objects like gauge-dependent correlators. In the 
non-Abelian case, we must therefore consider currents associated with additional global 
symmetries among the matter fields, such as the SU(2)r custodial symmetry of the SU(2)- 
Higgs theory. 
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j(t) < > j(0) 



FIG. 3. The leading-order perturbative diagram for (j(i)j(O)). 



One way to approximate correlators is by straightforward expansion in diagrammatic 
perturbation theory. Fig. ^| shows the leading contribution to a current-current correlator 
(j(t)j(O)). The result (for scalar fields) is 

V(m 3j (0)) = l5 ljq 2 [ vlln^e^ + il + n^e-^} 2 , (2.1) 

where u p is the energy of an excitation with momentum p, v p = V p c<jp is its group velocity 
(so v 2 = (p/ujp) 2 for Lorentz invariant continuum theories), and / is short for / d 3 p/(27r) 3 . 
(See appendix [A| for the case of fermionic charge carriers.) is the Bose distribution, 

^ = ( 2 - 2 ) 

where % has been shown explicitly for the first and last time in this paper. The case of 
classical field theory may be obtained by taking the h — > limit of the Bose distribution, so 
that 

1 + n w -»• n u -»• . (2.3) 

jjui 

In (|2.1|) , q 2 is the appropriate sum of squared charges, for the particular current under con- 
sideration, of different particle species, (q 2 is 1 for the number current of a single complex 
scalar. For non-Abelian currents such as (|1.2| ), q 2 should be multiplied by 2S a b.) For conve- 
nience, we will assume throughout that the currents of interest are normalized such that q 2 
is 0(1). 

We can split the result ( |2.1[ ) into an asymptotic constant q 2 C and a decaying transient 
q 2 D(t) by rewriting it as 



VimjM) = q 2 [C + D(t)\, (2.4a) 



with 

C 



f v 2 p 71^(1 + n Up ), (2.4b) 



D(t) = \ vl <e 2 ^* + (1 + n^fe- 2 ^ 1 . (2.4c) 



3 / 
'P 



Note that (j(t)j(O)) is real for classical field theory but not for quantum field theory, because 
in the latter case j(i) and j(0) do not commute. The imaginary part only appears in the 
transient q 2 D(t). 
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The asymptotic constant q 2 C has a simple physical interpretation: Eq. (|2.4b ) is precisely 



the result one would obtain for the expectation (j(O)j(O)) in an ideal Bose gas of particles 
(and anti-particles) rather than of fields. These particles may be regarded as classical (with 
definite position and momentum) except that they satisfy Bose or Fermi statistics. We 
shall refer to this as a kinetic theory description of the plasma. The current is gv for 
each particle. Because velocities of particles do not change with time in an ideal gas, 
(j(^)j(O)) = (j(O)j(O)) — <? 2 ( v (0) v (0)) is a constant for such a gas. Note that q 2 C is real, as 
it should be since is real and commuting in the kinetic theory approximation. 

The decaying transient q 2 D(t), on the other hand, can be seen to arise from an interfer- 
ence between particle and anti-particle contributions in field theory. Consider, for example, 
a complex scalar field <fi written in terms of creation and annihilation operators: 



[x 



I 4^ (V^ + 6i p e^ pi ) > (2-5) 

Jp y ZUl p 



where we have used non-relativistic normalization for a p and b p . A (space-averaged) number 
current in a continuum scalar theory would then be 

l(t) = y f x f V0 = J v p [JV P + N_ p + albl p e- 2 ^ + 6_ p a p e 2 ^*] , (2.6) 

where N p = a p a p and N p = b p b p are the number operators for particles and anti-particles, 
respectively (and J x is, of course, short for / d 3 x). Now consider the correlator (j(t)j(O)). 
The non-oscillatory terms in the current ( |2.6| ) correspond to the kinetic theory description 
and are responsible for the asymptotic constant discussed before. The oscillatory pieces 
of the current, involving interference between particle and anti-particle contributions, are 
responsible for the transient behavior of the correlator. These oscillations become irrelevant 
at large times because contributions from different momenta p decohere. 

If we treat the particles as massless, then the only frequency scale in the problem is T 
for the continuum theory, or the cutoff a~ l for the classical lattice theory. The time scale 
for the decay of the transient D(t) is therefore 1/T or a, respectively, and this is the origin 
of the transients at those scales in figs. |l] and ||. 

In the case ( |2.3p of classical field theory, for a current carried by bosons, the expres- 
sion ( |2.4| ) for the transient contribution directly gives -D(O) = C. Hence, in leading-order 
perturbation theory, the correlation decays by exactly a factor of 2 from its initial value, 
as shown in fig. 0. The same thing is true of the thermal contribution to the correlator in 



quantum field theory. The zero-temperature contribution to the correlator [n — > in (2.4)1, 
in the limit of negligible masses (so that u p = |p|), generates a 1/t 3 imaginary part for the 
correlator. 

The decay of the transient can be understood by analyzing the long-time behavior of 
D(t).[] Expressing (|2.4c|) in terms of a spectral weight, 



1 Long-time here means long compared to 1 /T (or a) but applies only to what can be seen in the perturbative 
approximation of fig. |j| It does not refer to the very long-time (t > 1/ '<7 4 Tln) physics shown in fig. [l| 
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D(t) = |Jf du; ^(w) [nl e 2wi + (1 + n^e" 2 * 

/oo 
dujp vv (uj)nle 2i "\ (2.7) 
-oo 

with 

p OT (u;) = / vl5(uj-ujp), (2.8) 



p 



converts the problem to a standard exercise analyzing the asymptotic behavior of a Fourier 
transform. (In Q2.71) , we used 1 + n w = — This asymptotic behavior is determined by 
non-analyticities of the integrand. If we ignore the mass m, then in the continuum case with 
uj p = p, the integrand is analytic in a strip about the real axis. Consequently, the long-time 
behavior decreases exponentially with an exponent determined by the imaginary part of the 
nearest singularity in the complex u plane (which is at ±2iiiT). So the decay shown in fig. [I] 
at t ~ 1/T is exponential. In fact, in this particular case, one can find an exact analytic 
result for all times, 

C = |T 3 , D(t) = [(MTt + iir) coth(2nTt) - l] . (2.9) 

For a classical lattice theory, the situation is somewhat different. In this case, ( |2.3|) and 
( P- 4 b| ) then imply that the steady-state value C is dominated by momenta of order the inverse 
lattice spacing 1/a rather than T, with the result that C ~ T 2 /a. Also, the density p vv (u) 
(2.8) has kink singularities generated from points at the edge of the Brillouin zone where 



v p = 0. These are the well-known Van Hove singularities which appear in the electronic 
density of states in crystals. Such kink singularities produce an oscillating power-law tail 
in Fourier transforms, with the frequency of oscillations determined by the location of the 
singularities. This is the source of the oscillations in the t < a transient shown in fig. |2]. The 
envelope of the oscillations decays like t _5//2 ; a detailed expression for the long-time behavior 
of D(t) is given in appendix |B[ 

Finally, let us go one step beyond the naive perturbation theory we have considered 
so far. It is well known that perturbation theory can be improved for high temperatures 
by resumming "hard thermal loops" f[l| . In particular, this means including the effective 
thermal mass in propagators. The thermal contribution to the mass m is O(gT) in the 



continuum, and 0(g^jT / a) for classical lattice theories. In the presence of a mass, there 
will be discontinuities in p vv (uj) at u — ±m. These give rise to an additional power-law 
contribution to the long-time behavior of D(t) that oscillates with frequency m. These 
oscillations have not been shown in figs. [3] and because they are sub-leading compared to 
the total correlation: the amplitude of the transient has already decayed to 0(g) by the 
time the mass oscillations begin. See the complete expressions in appendix [B] for details. 



III. KINETIC THEORY 

In the last section, we saw that the leading-order perturbative result for the correlator 
quickly settles down to the kinetic theory result for a non-interacting gas. In the kinetic 
theory picture, it is easy to see that the effect of interactions will become crucial at late times, 
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FIG. 4. A pictorial example of a typical fluctuation with net current j. 



and so simple (one loop) perturbation theory must break down.Q Currents are produced by 
fluctuations in the plasma. Imagine focusing attention on some region of space (with a size 
large compared to the inter-particle separation). Suppose that, at time zero, a few more 
positive charged particles happen to be moving left rather than right and/or a few more 
negative charged ones are moving right rather than left, as depicted in fig. [|, so that there 
is a net current j(0) to the left. As time progresses, the directions of the particles will be 
progressively randomized by collisions, causing the net current to decay and the correlation 
(j(t)j(O)) to approach zero. 

A quantitative analysis of the effect of interactions will be enormously simplified if one can 
indeed forget about quantum field theory for times large compared to 1/a or 1/T, and instead 
study the dynamics using kinetic theory. Kinetic theory applies to times and distances 
large compared to particle energies and momenta. In addition, the validity of kinetic theory 
requires that de Broglie wavelengths and scattering durations be small compared to the mean 
free path and mean free time between collisions, respectively. This ensures that particles 
can be viewed as propagating classically (with on-shell energies) between collisions which 
are independent, uncorrelated events. Whether these conditions are satisfied depends on 
the typical momenta and dominant scattering processes that contribute to a given physical 
problem. In the case at hand, it will be easiest to assume the validity of kinetic theory and 
return later to check the conditions. 



2 The same conclusion can, of course, be reached from a diagrammatic analysis, albeit with more effort. For 
long times, one finds that contributions comparable in size to the lowest order diagram can be generated 
from an infinite series of ladder-like diagrams which contain (near) on-shell singularities. Summing these 
diagrams yields the same result which is obtained from an effective kinetic theory description. For a detailed 
analysis, in the context of a scalar theory, see |12|. 
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FIG. 5. t-channel gauge boson exchange. The solid lines represent any sort of particles that couple to 
the gauge force. 

Mean free times 

To begin, it will be convenient to review some well-known basic scales associated with the 
times between collisions. Scattering among particles that carry gauge charges is dominated 
by t-channel exchange of a gauge boson, as depicted in fig. ||. Consider a collision between 
typical particles in the plasma, each with momentum of order T. For continuum relativistic 
particles, the differential cross-section is 

nda ~ng — 2 g T— (3.1) 

% u ' 

if we (temporarily) ignore plasma effects on the interaction. Here, £m = — Q 2 is the virtuality 
of the exchanged gauge boson, 9 is the angle by which one of the particles is scattered, and 
we've multiplied da by the density n ~ T 3 of scatterers. The time between single, large-angle 
{9 ~ 1) scatterings of particles is therefore 0(l/g A T). 

The total cross-section, however, is dominated by small-angle scatterings. Eq. ( |3.1| ) looks 
like it will lead to a linear infrared divergence in t, but plasma effects appear when y^M is 
comparable to the inverse Debye screening length ~ gT (corresponding to 9 ~ g). At this 
scale, Coulomb interactions are Debye screened, and magnetic interactions are suppressed 
by Landau damping of time-dependent, space-like magnetic fields.^ The damping of such 
magnetic fields is a reflection of the conductivity of the plasma and Lenz's law that a 
conducting system resists changes in the magnetic field.Q If all Coulomb and magnetic 
interactions became completely irrelevant for a/^m <C gT, then the mean free time for 
scattering through any angle would be order l/g 2 T and dominated by 9 ~ g. However, the 
plasma does not damp static magnetic fields, and so it does not completely damp nearly- 
static ones. Magnetic interactions remain significant in the small region of phase space in 
which the energy transfer (in the plasma rest frame) of the colliding particles is very small. 
This phase space suppression turns outQ to convert the linear small-tM divergence of ( |3.1|) 
into a logarithmic divergence dt M /t M when a/^m ~ gT. In non-Abelian gauge theories, this 



3 At the technical level, one incorporates these effects by including the one-loop thermal self-energy in the 
propagator of the exchanged gluon. 

4 We thank Guy Moore for pointing out this simple explanation. 

5 For details see, for example, ref. Jlg| l. 
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logarithmic divergence is then cut off by non-perturbative magnetic physics at a scale of 
\ftu ~ 9 2 T (0 ~ g 2 ). Hence, there is a logarithmic enhancement in the scattering rate and 
so a logarithmic suppression in the mean free time r sma n for any-angle scattering: 

Tsma11 ~ prlr? (3 ' 2) 

where "In", for non-Abelian theories, means ln[(gT) / (g 2 T)} ~ ln^^ 1 ). The mean free 
time for any-angle scattering of gluons in particular is typically referred to as the plasmon 
damping rate (in this case, for hard gluons). 

A different quantity of importance is the time it takes for the direction of a single particle 
to be changed by a large angle (0 ~ 1). As discussed above, this can happen by a single, 
large-angle scattering in a time of order \j g A T . However, it can alternatively happen by a 
succession of scatterings through a small angle AO. The time between such scatterings is of 
order (A0) 2 /g 4 T, provided that y/tu ~ gT (or A9 ^ g) so that there are no suppressions 
by the screening and damping effects discussed above. Treating successive scattering as 
a random walk, order \/{A0) 2 such scatterings are needed to build up a total deviation 
of ~ 1. So a succession of scatterings by AO also takes a time of order l/g 4 T to give 
large-angle scattering. Since all possible decades of AO contribute equally, the net result is 
a logarithmic enhancement in the large-angle scattering rate, r^g ~ g 4 T J d(A0)/(A0). The 
mean free time Tu g for large-angle scattering is then 

Tbis ~ ^ln (3 ' 3) 

This time, the logarithm is cut off by T in the ultraviolet and gT in the infrared, but it's 
still ln(g _1 ). The non-perturbative physics associated with the scale g 2 T in non-Abelian 
plasmas is not involved at leading order. 

For classical lattice theories, the above analysis can be repeated, but the typical mo- 
menta of particles is a -1 rather than T, the density is n ~ Ta~ 2 , the Debye mass scale 
is (3^ l ^a~ l , and the scale for non-perturbative magnetic physics in non-Abelian gauge 
theory is P^a" 1 , where the lattice coupling P^ 1 is order g 2 Ta in physical units. This 
results in Tbig ~ P^ 2 a/ \n(P^) ~ [g 4 T 2 a ln^)] -1 for the large-angle scattering time and 
Tsmaii ~ Pj^a/hx ~ ((? 2 T 2 aln) _1 for the any-angle scattering time, where "In" is ln(/?L) in 
non-Abelian gauge theory. 



The correlator 

So which scattering time controls the decay of the (j(t)j(O)) correlator? Return to the 
picture of a typical fluctuation at time zero that has a net current j(0), illustrated in fig. |j. 
The total current j = Yl s 3^ is the sum of currents = g^v^ of the individual particles 
(labeled by s). Any one of those individual currents will only change substantially when 
the velocity v^ s ' changes substantially^ That means the relevant time scale for the random- 



6 Recall that we have assumed that the current j is gauge-invariant, and so the corresponding symmetry 
commutes with the gauge interactions. As a result, the are not changed by small-angle scatterings due 
to gauge-boson exchange. 
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ization of the current is the mean free time Tbi g for large-angle scattering, rather than the 
mean free time T sma u for any-angle scattering. This is the physics responsible for the period 
of exponential decay starting at time l/(g 4 Tln) in fig. [I] [and time l/(g 4 T 2 aln) in fig. 0. 

There has been confusion in the literature on the above point. In ref. |[10|| , for example, 
it was suggested that the way to understand the decay of the correlators measured by Tang 
and Smit is simply to include thermal widths (the imaginary part of the self-energy) in the 
internal propagators of the leading-order perturbative diagram, fig. |3|. The imaginary part 
of the self-energy measures the rate for a particle to scatter into any other state and so 
includes very small-angle scatterings — that is, the width is of order r^ n . The inclusion of 
widths in the diagram of fig. |3| makes the internal propagators decay on a time scale of r sma u 
and therefore appears to predict that the correlation decays on a time scale of r sma n instead 
of Tbig. The problem with this reasoning becomes apparent if one realizes that it could be 
applied to any correlator at all, including the total charge correlator (j°{t)j°{0)) (where 
j°(t) = J x j°(x, t)). By conservation of charge, the latter does not decay at all with time. 
The problem is that the decay seemingly generated by resumming self-energies into fig. [3] is 
canceled by other, higher-order diagrams in the resummed perturbation theory. An explicit 
example of such cancelations for QED plasmas can be found in ref. [|14| . This highlights the 
dangers of using resummed perturbation theory for intuition about physics at large times. A 
(highly non-trivial) discussion of how kinetic theory may be seen to arise from perturbative 
diagrams in scalar theory may be found in ref. Jl2| . 

We shall now outline how the rate of exponential decay (for times of order r^g) can be 
calculated quantitatively from kinetic theory. One starts with the Boltzmann equation for 
the probability distributions /(x, p;t) of particles (and anti-particles) in phase space as a 
function of time. The Boltzmann equation has the form 

(d t + v p -V x )f = C[f]. (3.4) 

On the left is a convective time derivative of the distribution function. The right hand 
side is a collision integral representing the net rate at which particles are scattered into a 
momentum state p minus the rate at which they are scattered out. For 2 — > 2 scattering, it 
has the form 

C[f] = I |M pp ^ kk f [/ p /p,(l + / k )(l + fv) - hMl + / P )(l + /p')] , (3.5) 

Jp'kk' 

where M. is the scattering amplitude (and indices labeling particle species are suppressed). 

We are interested in the response in time of a fluctuation away from equilibrium in the 
(space-averaged) current 

m = V-' I I 9v p /(p,x;t). (3.6) 

Jx Jp 

The response to fluctuations can be analyzed by linearizing / about the equilibrium distri- 
bution n(u; p ).[j The resulting linearized Boltzmann equation has the form 



7 A classic discussion of most of these techniques can be found in chapter IV of ref. [[l5| . See also ref. Jl6| 
for a somewhat related application of calculating the shear viscosity of hot QCD. 
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(<9* + v p -V x )5f «C[n] Sf, 



(3.7) 



where C[n] is a linear integral operator acting on Sf. The decay of the space-averaged 
component of Sf is then determined by the eigenvalues of C[n}. The eigenvector with smallest 
eigenvalue describes the most slowly relaxing fluctuation (with a decay constant given by the 
eigenvalue), and so dominates the long-time behavior in this approximation. The operator 
C[n] has zero modes corresponding to fluctuations in total charge, energy, and momentum, 
but in the symmetry channel corresponding to the total spatial current (C odd, CP even 
fluctuations) the smallest eigenvalue of C[n] is non-zero. In our companion paper M, we 
show how to calculate this decay exponent explicitly in leading-log approximation. The 
result is indeed 0(r^), as claimed earlier. 

Having realized that the relevant time scale is 0(Tu g ) rather than 0(r sma n), we can now 
check the validity of kinetic theory. We need to verify that the de Broglie wavelengths 
and scattering durations are small compared to the time between these scatterings. The 
fluctuations in the current are dominated by hard particles with momentum T and de Broglie 
wavelength 1/T. The dominant processes that contribute to the large- angle scattering time 
O(rbig) are, as discussed earlier, scattering events with momentum transfer between gT and 
T. The time between such individual scatterings is l/g 2 T and l/g 4 T, respectively, and so 
we need kinetic theory to be valid for time scales as small as At ~ l/g 2 T. The scattering 
duration for scattering with momentum transfer is ordeif] l/\q — Qo\ ~ \j \[Q 1 % 1/gT. 
This, and the de Broglie wavelength, are indeed small compared to At for weakly coupled 
theories, and so kinetic theory is appropriate.^ 



Large Nf scaling 

We digress momentarily to mention the scaling of the characteristic time Tbig for the decay 
of current correlations when there are a large number Nf of equivalent charged fields. As 
usual, the gauge coupling is assumed to scale as l/Nf so that A = g 2 Nf is fixed as Nf — > oo, 
and the inverse Debye screening length ~ a/AT has a finite large Nf limit. Repeating the 
earlier discussion, one finds that the mean free time for large-angle scattering grows linearly 
with Nf, 

Nf 

Tbig ~ A*rin(A-i) ■ (3 ' 8) 

This means that there is a non-uniformity between the large time and large Nf limits, 
with limAr ^ oo lim t _ +oo (j(t)j(0)) = 0, while ^m t - too ]im N .- ¥OO (j(t)j(0)) ^ 0. Consequently, 



8 This estimate holds for typical scattering processes for which q — Qq and q + Qq are of comparable 
magnitude. For highly collinear scattering events, \q — Qq\ <C \q + Qo\, and the scattering time can be much 
larger than 1/gT. However, the phase space of these exceptional scattering events is sufficiently restricted 
that they make a negligible contribution to Tbi g (or any of the physics under discussion) . 

9 If the relevant time scale for the physics of interest had been 0(r sma ii) rather than 0(Tu g ), then kinetic 
theory would be inappropriate (except at leading log order): it would break down for momentum transfers 
of order g 2 T that contribute importantly to 0(T sma n). 
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calculations which use 1/Nf expansion techniques |3j]], valid for arbitrary fixed times t, will 
not be able to see the eventual decay of correlations occurring on time scales which grow as 
0(JV». 



IV. LONG-TIME HYDRODYNAMIC TAILS 

In the last section, we argued that the velocities, and therefore the currents, of individual 
particles get randomized on a time scale of Tu g - A picture of this randomization, as a result 
of successive large-angle scatterings, is shown in fig. |]a. Based on it, one might expect 
the current correlations to show exponential decay at arbitrarily long times (exactly as 
predicted by the linearized Boltzmann equation). This prediction, however, is wrong, as 
was first discovered in 1970 in numerical simulations of a gas of classical billiard balls |TJ| . 



It was found that if one follows the velocity of a particular billiard ball s, the velocity 
auto-correlation function (v"(t)v^(0)) has power-law decay at large times. Such power- 
law decay is a consequence of the presence of hydrodynamic fluctuations with arbitrarily 
long relaxation times. Q Fig. |6|b illustrates the relevant physics. Imagine a fluctuation in 
which there is a very long and wide column of fluid that is in local equilibrium, but has 
a net average velocity v nct with respect to the rest of the fluid. After several large-angle 
scatterings, the velocity auto-correlation (v' s '(t)vW(0)) of particles well inside this column 
will not fall to zero, but will approach the non-zero local equilibrium value (v^ s ^) 2 = v 2 et . 
This remaining correlation will only decay on a time scale determined by the diffusion of the 
particles out of the column of moving fluid. (This includes two related effects: the particle s 
under consideration could leave the column, or the column's velocity could spread out and 
degrade by the transverse diffusion of the other particles. The latter process is parameterized 
by the shear viscosity of the fluid.) Because diffusion is a random- walk process, the time it 
takes for the remnant velocity correlation v 2 et to decay is proportional to the square of the 
transverse size L of the column: t ~ (L/v) 2 /tu s , where v is the typical particle speed. 

Now that we've sketched the basic point, it is no longer necessary to consider fluctuations 
of infinite spatial extent like the fluid column in fig. [|b. In general, after an arbitrarily long 
time t 3> Tbig, the velocity auto-correlation of a given particle s will be reduced to roughly the 
square of the average velocity of the particles within causal reach of particle s by diffusion — 
that is, the particles within a volume a radius L ~ v^/r^ g t. The number of such particles 
is N ~ nL 3 , and this average velocity is then order v/y/N. So, at late times, 

(vW(t) ■ v ( s >(0)> ~ ^ ~ [nv (r big t) 3 / 2 ]^ (t » r big ). (4.1) 

This t~ 3 / 2 power-law decay is precisely what was measured in simulations of billiard-ball 
gases .0 

Now turn to our case of interest: the decay of the current- current correlation (j(t)j(O)). A 
long-wavelength fluctuation only in the velocity, as shown in fig. |6|b, is not enough to create 



10 For a short discussion in the early literature, see ref. For a review of the entire subject, see ref. [^0| 

11 In d spatial dimensions, this becomes a t~ d / 2 power law decay. 
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FIG. 6. (a) The randomization of a particle's velocity by successive scatterings in the plasma, (b) The 
same process in the background of a large column of fluid moving with net velocity. In this case, the particle's 
velocity is not completely randomized but is biased by the net velocity v nct . 

a long-lived net current. We must additionally have a roughly coincident long- wavelength 
fluctuation in the net charge, as shown in fig. [7[ This is a long-lived fluctuation in the 
current because, due to charge and momentum conservation, neither the net local charge 
nor the net local velocity can dissipate except by diffusion, and diffusive dissipation takes 
arbitrarily long if the wavelength of the fluctuation is arbitrarily long. (This is in sharp 
contrast with the situation shown in fig. where there is a net current but no net local 
velocity or charge.) Because the current in fig. [7] is the product of two small fluctuations 
from equilibrium — the charge fluctuation and the velocity fluctuation — it is a second-order 
effect in fluctuations and will not be seen in a strictly linearized analysis of the problem. 
This is why the linearized Boltzmann equation discussed in section [If| does not yield a 
very-long-time power-law tail.0 

We will show in a moment how to use the equations of hydrodynamics to see that the 
current correlator decays as t~ 3 ^ 2 , but first we give a very rough and heuristic picture of 
this result. Let j ~ qv be the typical absolute magnitude of the current of an individual 
particle. Then the total current at time zero is of order jy/nV, where nV is the total number 
of particles in all of space. This correctly reproduces V(j(0)j(0)) ~ nq 2 v 2 ~ T 3 , where j is 
the space-averaged current. Now consider the (j(i)j(O)) correlation for times large compared 
to Tbig. In the case of the billiard-gas velocity correlation, all that was important was the 
average velocity of a cell roughly of size L ~ Vy/r^g t. For the current correlation, what's 
important is both the average velocity and net charge of such a cell. These are of order 
v/y/nV* and qy/nLP, respectively, where nL 3 is the number of particles in the cell. This 
gives a contribution to the current- current correlator of order qv for each cell. Since there 
are iV cells ~ V/L 3 such cells, each with randomly directed contributions to the correlation, 



12 It would be interesting to analyze in detail the Boltzmann equation treated to second-order in fluc- 
tuations, as it bears on our problem, and verify explicitly the time scale (discussed below) at which the 
linearized Boltzmann equation must break down, but we have not done so. 
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FIG. 7. Coincident hydrodynamic fluctuations in both the velocity and charge density of the fluid, 
giving rise to a long-lived fluctuation in the net current. (Two columns are shown, with parallel currents 
but opposite momentum and charge densities, merely to emphasize that large scale current fluctuations may 
occur with no change in total charge or momentum.) 



the total will be of order 



N cells 2 q 2 /v 



1 



(t > r big 



(4.2) 



We can get a rough idea of when this power-law decay takes over from the exponential 
decay discussed in the last section by equating the correlator in the two approximations: 



y^g — icross/T big 



1 



(^big Across ) 3 ^ 

implies 

Across ~ ibig hi(r big T) ~ l/(g 4 T). 
At this point, the amplitude of the correlation is 

vw™ ss ) -'M) ~ {g A T?^'\g- 1 ), 

as indicated in fig. |l|. 



(4.3) 



(4.4) 



(4.5) 



Hydrodynamic equations 

To be a little more systematic about the analysis of long-time tails, we will now repeat 
our discussion within the context of the natural effective theory for describing long-time 
excitations in plasmas: hydrodynamics. In general, hydrodynamics is the description of long 
wavelength fluctuations of conserved densities, i.e., energy, momentum and charge densities. 
Long wavelength fluctuations in these quantities can persist an arbitrarily long time because 
they can relax only by diffusive processes which transport the conserved quantities over long 
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distances. Consider a simple system with one conserved charge q in addition to the energy 
and momentum. Let e = T 00 , 7r, = T l °, and p = f be the energy, momentum, and g-charge 
densities. Hydrodynamics consists of conservation laws, 

d t e = -V • 7T, (4.6a) 
dtTTi^-VjTij, (4.6b) 
«9 t p = -V-j, (4.6c) 

and constitutive relations which express the spatial currents in terms of the conserved densi- 
ties in the hydrodynamic limit. In this example, we need constitutive relations for the stress 
tensor and the spatial charge current j. The constitutive relations may be parameter- 
ized by writing down, in a small-momentum expansion, all possible terms consistent with 
symmetries (including boost invariance) of the resulting equations. Typically, one expands 
only to first-order in momentum. For our purposes, it will be least cumbersome to also 
consider an expansion in fluctuations 5e, n, and p away from global equilibrium.^ Then, 
for example, terms allowed by symmetry for linearized hydrodynamics are 

Tij = 5ij (a + b5e- 7 ? V • 7r) - ~/ v (V^- + V,-^ - §<%V ■ 7r) , (4.7a) 
j = -D q Vp, (4.7b) 

to first order in gradients. The constants are usually expressed in more physical terms as 
a = P, b = vl , j v = r]/ (e + P), and 7^ = (/ (e + P), where P is the pressure, v B is the speed 
of sound, 7] and ( are the shear and bulk viscosity, and (e + P) is the equilibrium enthalpy 
density. D q is the diffusion constant for the conserved charge q. 

In the linearized approximation, the space-averaged current is automatically zero: 

j(t)«- / A,Vp(x,t) = 0. (4.8) 

To study correlations of the space- averaged current, we will therefore need to consider the 



next higher order terms in the current constitutive relations ( f4.7|) . Incorporating the next 
term consistent with symmetries^ produces, 

j = — D q Vp + cp 7T + (higher order). (4.9) 

In fact, the constant c is determined by boost invariance, and this expansion is just 

j = — D q Vp + p v + (higher-order), (4-10) 

where 

v = + (higher-order) (4-H) 

(e + P) 



13 As everywhere in this paper, we are assuming that there is no chemical potential(s), hence the equilibrium 
change density is zero. We also choose to work in the rest frame of the fluid, so that iz vanishes in equilibrium. 



14 Including charge conjugation, under which 



'ti- 



ll 



is the local fluid velocity. This has been a rather roundabout way to realize that there is 
a contribution of the form pv to the current, but this methodical approach will be useful 
when we move on to discuss long-time tails in lattice theories. 

We can now analyze the hydrodynamic decay of current correlations by using the second- 
order constitutive relation ( |4. 10| ) for j to express 

(j(f).j(0))~(/ p(x,t)v(x,t). / p(y,0)v(y,0)> 

= ( j p(k, t)v(-k, t) ■ j p(-k', i)v(k', 0)), (4.12) 

where p(k, t) denotes the spatial Fourier transform for p(x, t), etc. One may now use lin- 
earized hydrodynamics to compute the correlation (which factorizes at this order), 

-J(0)> ^ / (p(k,t)p(-k',0))(v(-k,t)-v(k',0)). (4.13) 

Solutions to the hydrodynamic equations tell one about the decay of correlations. The 
linearized p equation, dp/dt = D q V 2 p, gives 

p(k,t) ~ p(k,0)e- D « k2t , (4.14) 

which should be interpreted in the present context as meaning 

(p(k, t)p(-k', 0)> ~ (p(k, 0)p(-k', 0)) e~ D « kH . (4.15) 

To find the equal-time fluctuation spectrum (p(k, 0)p(— k', 0)), we have to step back from 
hydrodynamics and realize that in position space the charge density p is only correlated over 
very small distances (of order 1/T). As far as the long distance scales of hydrodynamics are 
concerned, this correlation can be approximated as local: 

(p(x,0)p(y,0))~ X ^ (3) (x-y), (4.16) 
where \q is the charge susceptibility, 

Xq = j (p(x,0)p(0,0)>. (4.17) 

Consequently, 

(p(k,t)p(-k',0)> ^ Xq 6^(k-k')e- D ^ 2t . (4.18) 

The corresponding solutions for the behavior of velocity v(k, t) (or momentum density) fluc- 
tuations split into the two cases of transverse (v _L k) and longitudinal (v || k) perturbations. 
The solutions in the two cases behave as 

v ± (k, t) ~ v ± (k, 0) e~ 7 " fc2 *, (4.19a) 

v,| (k, t) ~ v,| (k, 0) cos(v s kt)e-^ sk2t/2 - -r^y 8e(k, 0) sin^t) e" 7sfc2t/2 (4. 19b) 
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Transverse fluctuations ( |4.19a| ) relax diffusively at a rate controlled by the shear viscosity, 
while longitudinal fluctuations (|4.19b| ) are propagating sound waves which decay at a rate 
governed by the sound attenuation constant 7 S = + 7^. The net result for the current- 
current correlator ( |4.13| ) is 



where the susceptibility Xv is 



2 / e -iD q+ ~ h )kH + j CQ ^ vM) e -(D q+ls /2)k'h 



k 



(4.20) 



X" = 3(e + P) 2 JJ*(* °) • ""(0. °))- ( 4 - 21 ) 

The first integral gives the t~ 3 / 2 behavior discussed previously. The second integral decays 
exponentially in time due to decoherence in the spectrum of sound waves. The long-time 
result is then 



[47r(D q + lv )t} 



V(i(t) ■ j(0)> * - - — - • (4.22) 



The susceptibilities Xq an d Xv are equal-time quantities and are easily computed using 
perturbation theory in the original quantum field theory. They are order T 3 and T -3 re- 
spectively. The constants D q and 7^ are related to diffusion and are 0(tu s )- They can be 
evaluated precisely in any given theory from the linearized Boltzmann equation. Techniques 
for calculating such coefficients, under the simplifying assumption that ln(g _1 ) ^> 1, can 



be found in ref. [16], which computes the shear viscosity 77 (and so 7^) for QCD in this 



approximation. In any case, the order of magnitude of the correlator at long times is 

,/(j(t) ' j(0)> ~(^' (423) 

in agreement with our previous estimate 



Not-quite-hydrodynamics on the lattice 

Now let's try the same analysis for theories on a spatial lattice. Because the lattice 
breaks continuous translation invariance, there is one crucial difference between the lattice 
and continuum theories: momentum is not an extensive, additively conserved quantity in 
the lattice theory. More specifically, momentum is only conserved modulo 2ir/a on a (simple 
cubic) lattice. In the continuum, boosting an equilibrium state yields another equilibrium 
state: A uniformly boosted plasma, with non-zero momentum density 7r and total momen- 
tum irV, is a stable equilibrium state. On the lattice, however, the total momentum irV 
will decay away due to microscopic 2^2 particle collisions^ that violate momentum by 



15 Depending on the details of the lattice Hamiltonian, 2 — ► 1 and 1 — > 2 processes may also be possible, 
in which case the time scale discussed below will be even shorter. For the canonical Kogut-Susskind Hamil- 
tonian, however, 2 — > 1 and 1 — > 2 processes are forbidden by energy and lattice momentum conservation. 
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multiples of 2-7r/a in one or more lattice directions. (These are just the umklapp processes of 
solid state physics.) This 2 — > 2 process is dominated by momenta, and momentum trans- 
fers, of 0(1 /a). The corresponding time scale is the same as that for individual large- angle 
scatterings: l/0(g A T 2 a).f^\ As a result, large-scale fluctuations in momentum may decay 
on the same time scale as other microscopic scattering processes; consequently, momentum 
density should no longer be treated as one of the fundamental independent variables for 
"hydrodynamics" on the lattice. 

So what does the appropriate long-distance and long-time theory look like? The conser- 
vation laws ( f4.6p are now reduced to 

d t e = -V-7r, (4.24a) 
d tP = -V-j, (4.24b) 

and we now need an effective constituent equation for the momentum density 7r, as well as 
j, in terms of just energy and charge densities. The lowest-order possibility for 7r is 

7r = —D e Ve + (higher order). (4-25) 

The lowest-order possibility for the current j(x, t) that will give a non- vanishing contribution 
to the spatial average of j is 

j = —D q Vp + Ci 5eV p + c 2 [p, Vp] + (higher order). (4.26) 

The commutator term [p, Vp] appears only if the current j is associated with a non-Abelian 
global symmetry (as is the case for Tang and Smit's W correlator). 

Following the procedure we used for the continuum case, the long-time tail of the current 
correlator is then given by 



v(j(t) • j(o)> * c\ XeXq [ e e -^+ D ^ + 4xi j fcV 

[A7T{D £ + D q )t] 5 / 2 [SttIVF 2 ' 



2D„k 2 t 



(4.27) 



where the susceptibility Xe is defined analogously to Xq (|4.17|) , with p — > 5e, and is better 
known as T 2 times the specific heat (more precisely, T 2 times the heat capacity per unit 
volume). 

As in the continuum quantum theory, the susceptibilities can be computed in perturba- 
tion theory, and the diffusion coefficients could be calculated from the linearized Boltzmann 
equation. The susceptibilities are of order 



X 



~ TV 1 , X e ~ T 2 a" 3 , (4.28) 



16 It would be convenient to work in lattice units as then everything can be expressed in terms of the single 
variable /3l ~ (g 2 ^^ 1 instead of g, a, and T separately. In lattice units, the current should be scaled with 
an extra power of g 2 , so that Jiat(tlat) ~ <7 2 a _3 j(aiiat)- However, we will stick to physical units to facilitate 
comparison with previous continuum results. In such comparisons, one should think of 1/a as being order 
T, since these are the UV momentum scales in the two cases. 
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and the diffusion constants are 



D q r big ~ (g 4 T 2 a In)" 1 , A ~ ((/^o)" 1 . (4.29) 

D e does not have a log suppression because the diffusion of energy density is controlled by 
umklapp processes which, as asserted earlier, are dominated by large-momentum transfer 
processes.^ (Log suppressions in other diffusion coefficients came, in contrast, from small- 
momentum transfers — that is, from small- angle scattering.) 

The final ingredients are the coefficients c\ and C2 of the constitutive relation ( [4.26 ). 
Physically, these characterize how fluctuations in j correlate with fluctuations in e and p. 
For simplicity, let's assume the symmetry associated with j is Abelian and ignore c^. We 
can then compute c\ from a ratio of equal-time correlators, 

k-o k 2 (|e(k,0)p(-k,0)| 2 > 1 ; 

1 ' hm^-(j(0)6(k,0)p(-k,0)), (4.31) 



XqXe k^Ok 2 

as may be verified by substituting ( |4.26| ) into Q4.3CQ . The equal-time three-point correlator 
can be computed, with resummed perturbation theory in the original field theory, from the 
diagram shown in fig. |8|. This diagram is dominated by momenta of order the thermal mass 
m ~ g^jT/a and yieldsQ 

Cl ~^ln(— ). (4.32) 
1 \maj 



The order of the long-time tail ( 4.27| ) is then 



The cross-over between exponential decay and the hydrodynamic tail occurs at 

T 2 a 4 ln 2M\ 

T 2 a e~ <cross / Tbig (4 34) 

[(A + Across] 5 / 2 ' 1 J 

analogous to ([4.3|) . This time is 



17 O ne wa y to understand how umklapp processes determine D e is to return to the continuum case of eqs. 
( 4.6b| ) and (4.7a). Schematically, tt ~ — w^Ve+7VV7r. Now add the effect of momentum decay on the lattice 

umk^ 



with a time scale i U mk to get, schematically, tt ~ —v^Ve + 'jW-k — t uTV. In the hydrodynamic limit of 



arbitrarily small frequencies and momentum, the ~k and VV-7T terms can be ignored, giving 7r ~ — WgtumkVe. 
Hence D e ~ v^t umk . 

18 When working in physical units, an easy way to get the order of a diagram in the classical theory is to 
include a power of T for each propagator and then make up dimensions with m if the diagram is infrared 
dominated or a -1 if UV dominated. 
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FIG. 8. Leading-order diagram for computing the (jep) correlation. 
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at which point the amplitude is 

Vm r ossH(0)>~2 20 T 12 a 9 In 2 



(Abelian symmetry) 



/3 L 12 ln 2 (AJ. 



(4.35) 

(4.36a) 
(4.36b) 



The scaling of this amplitude like the tenth power of /3£ 1 relative to the t = correlation 
(ignoring the logs) means that the hydrodynamic tail does not appear on a fine lattice 
(/?l ^ 1) until the correlation is very small. 

For a non- Abelian j, we need to consider the contributions to the hydrodynamic tail from 
C2 as well. These may be worked out similarly to the c\ contributions above and also give a 
t -5//2 tail. However, the C2 contributions give a much larger amplitude for the tail. One can 
easily find that 



c 2 



_i V 2 k-(trj(0)[p(k,0),p(-k,0)]) 

2 k^o k2(tr[p(k,0),p(-k,0)] 2 ) 
a 2 



mT ' 



and 



V(Ktcro SS )-m)~9 18 T n a 8 In 5 / 2 



(non-Abelian symmetry) 



(4.37a) 
(4.37b) 

(4.38a) 
(4.38b) 



as indicated in fig. 0. 



V. FINITE VOLUME 

So far, we have only considered the behavior of the correlator in infinite volume. Real 
lattice simulations are, of course, done in finite volume. How does this change things? 
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g*T 2 a In 
time > 

FIG. 10. Same as fig. |2|, but in a small volume L 3 (L < (g 4 T 2 a) _1 ). In here means ln(/3L) if L > m 1 

and ln(L/a) if L < m _1 , where m is the Debye mass m ~ gT l t 2 a~ 1 / 2 ~ /3 L ^a -1 . (It is assumed that L is 
still sufficiently large so that the Poincare recurrence time is far greater than the time scales of interest.) 
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First, consider the case where the volume is large compared to the scale Tbi g for hydro- 
dynamics. In this case, the discussion of the last section still applies, but now the largest 
fluctuation is limited to a size L, where L is the linear extent of the lattice. The largest time 
scale for dissipation is then ti ~ L 2 /ru g - After this time, the correlation will again decay 
exponentially. This can be seen more explicitly from fl4.27| ). The integral over k is now a 



sum, and the smallest k is 2n/L. The very-very-large time behavior of the propagator is 
then 

V(m -j(0)> - ± \c\ X ,X q e-^ +D ^ /L)H + 4x 2 q e-* D ^ H 



(t > L7r big ) (5.1) 



for an L 3 simple cubic lattice. The resulting behavior for the correlator is depicted qualita- 
tively in fig. ||. 

On the other hand, if the lattice is smaller than r^g, then the hydrodynamic power-law 
behavior will never develop. In this case, the correlator behaves as shown in fig. [Hj 



VI. WHAT ARE CLASSICAL SIMULATIONS GOOD FOR? 

Real-time simulations of classical, thermal, lattice gauge field theory are practical — one 
merely needs to integrate the classical equations of motion forward in time for a sample 
of Boltzmann-distributed initial conditions. Unfortunately, no comparable practical, non- 
perturbative method for studying real-time dynamics in thermal quantum field theories is 
known. For the specific purpose of computing the topological transition rate in hot non- 
Abelian theories, it was argued in ref. that the relevant degrees of freedom are essentially 
classical, and therefore it should be possible to determine the transition rate in the quantum 
theory by studying a hot classical lattice theory. In @ it was pointed out that the topological 
transition rate is sensitive to the damping rate for gauge field fluctuations with 0(g 2 T) 
momenta, and that this damping rate depends on scattering off the thermally-distributed 
hard excitations. In the quantum theory, the Bose distribution cuts off the distribution 
of hard excitations at momenta of 0(T), while a classical lattice theory instead has an 
0(1 /a) cut-off from the finite lattice spacing. Consequently, the damping rates for soft 
0(g 2 T) modes differ between the quantum and classical theories. Nevertheless, because 
these damping rates are perturbatively calculable [[/J, one can still, in principle, determine 
the quantum transition rate by applying a calculable correction factor to the (lattice spacing 
dependent) result of a classical lattice theory -M 



One might have hoped that thermal classical lattice theories, and hot quantum theo- 
ries, would generate the same behavior (up to calculable correction factors) for many other 
observables probing low-frequency, real-time response. The example of the current-current 
correlator studied in this paper shows that this cannot be universally valid. At short times 



19 We are glossing over the difficulties caused by the loss of rotational invariance in the classical lattice 
theory and resulting anisotropy in the classical damping rates However, this does not change the fact 

that essentially equivalent physics controls the dynamics of the relevant p ~ g 2 T, u> ~ g 4 T fluctuations in 
the quantum and classical lattice theories. 
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(t ~ 1/T or t ~ 1/a, respectively) the very different behavior of the current correlation 
in the quantum and classical theories (illustrated in figs. |l] and ||) is surely no surprise; at 
these times the correlation is sensitive to the precise form of the cutoff in the distribution 
of hard excitations, which is completely different in the two theories. However, the decay 
rate which characterizes the exponential regime is also sensitive to scattering rates off hard 
excitations, and no simple quantitative relation connects the resulting values for the classical 
and quantum theories. Finally, even the exponents of the long-time power-law tails differ 
due to the very different nature of hydrodynamic fluctuations in the quantum (continuum) 
and classical (lattice) theories. 

Naturally, the same issues arise in the analysis of virtually any other correlation function, 
although the details will differ. For example, the correlator of the gauge-invariant scalar 
operator / x $t(x)<I>(x), also studied by Tang and Smit 0, has much smaller short-time 
(t ~ 1/a) oscillations in the lattice theory because the one- loop perturbative contribution to 
this correlator is IR dominated, in contrast to the UV dominance of the current correlator. As 
discussed by Tang and Smit, the most noticeable feature is a decaying oscillation associated 
instead with the thermal mass, at time scales t ~ 1/m ~ l/gT.0 Since has a non- 
vanishing equilibrium expectation value, the two-point correlator of this operator does not 
decay to zero, but will instead approach the square of the (perturbatively computable) 
expectation value. Just like the current-current correlator, the time dependence of the 
correlator will show a plateau in the free kinetic regime followed (in infinite volume) 
by initially exponential and later power-law relaxation to the non-zero equilibrium value. 
And, just as was shown for the current correlator, quantitative aspects of this relaxation 
will differ between the continuum quantum theory and the corresponding classical thermal 
lattice theory. 

Some readers acquainted with the role of classical lattice simulations in studying the 
rate of high-temperature baryon number (B) violation may be uneasy. The time scale of 
B- violating processes is argued to be t ~ l/g^T at high temperature [fj]. But that is the very 
time scale where, for the case of current-current correlations, we've argued lattice simulations 
begin to fail to even qualitatively reproduce the continuum answer. There is, however, a very 
important difference between these two cases. At t ~ l/g 4 T, the correlator of the space- 
averaged current is dominated by collective fluctuations of hard modes (p ~ T), where the 
wavelength L of the collective fluctuation is order l/g 4 T. That is, large times correspond to 
equally large distances, and one is on the edge of the hydrodynamic regime (t, L 3> (g 4 T) _1 ), 
which the lattice fails to reproduce. B violating processes, in contrast, are dominated by soft 
modes (pCT). The soft modes of interest have momenta p ~ g 2 T, and the fluctuations of 
interest have wavelength L ~ l/g 2 T. It is these soft modes which should have an effective 
classical description, and the physics is different than the hard-mode physics that dominates 
hydrodynamic behavior. Nonetheless, one might still fear trouble^ because, as argued in 



20 Actually, the scaling of time scales is somewhat complicated in Tang and Smit's simulations, because 
they not only vary /?l ~ {g 2 Ta)~ l but vary the zero-temperature Higgs mass as well. In this paper, we have 
routinely treated the zero-temperature Higgs mass as negligible. 

21 That is, trouble more disastrous than what could be accommodated by some sort of perturbative 
matching. 
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ref. 0, the dynamics of the soft modes is affected by hard modes with which they interact. 
Does the lattice do anything terrible to the dynamics of these hard particles, and hence 
indirectly to the dynamics of the soft modes, on a time scale of l/g 4 T? This question was 
addressed by Huet and Son in ref. ||. The answer is that it doesn't matter. The time for 
a hard particle to fly across the region of a B violating process is given by the size of that 
region: L ~ l/g 2 T. What happens subsequently to the hard particle is irrelevant. Since 
l/g 2 T is small compared to Tbi g , the hard-particle collisions which produce hydrodynamic 
behavior are not relevant to the process of B violation. 
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APPENDIX A: FERMIONIC CHARGE CARRIERS 

If the global charge associated with the current j is carried by fermion fields, then the 
one-loop diagram fig. R] gives 



P 4o;2 



| 5ij q 2 



tr 



T< ((! A 
((1 -nu p ) AZ p e 



+ e -wt _ n A e 



V' 



iujt — n A + p iujt 



(Al) 



^ P ) + <e 2 ^ + (1 - n^fe- 2 ^ 
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P 



n u p e p + (1 - n Up ) e 



(A2) 



where^ Ap = ~^ij> + m and p° = 



\j[eP^ + 1) is the Fermi distribution. As in 



the bosonic case, the prefactor q 2 represents the sum of squared charges and, because of the 
two spin states, our convention is q 2 = 2 for the number current of a single Dirac fermion. 
The second term of the result (|A2|) gives a subleading contribution (suppressed by m/T) 
and will be neglected. Separating the initial transient from the asymptotic constant, as in 
(P]4k), leads to 



C =1 v l ~ n ^ P )' 



D(t) 



2 2iu p t 



+ (1 - n^fe 



(A3a) 
(A3b) 



At t — 0, one may easily see that Re -D(O) = — 2C. Hence, the current-current correlation 
rises from — C, crosses zero, and levels out at C. For times large compared to 1/T, the 



22 For those who prefer the H metric convention over 

be normalized so that Ap = ±]/> + m. 



the 7 matrices would conventionally 
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qualitative behavior is identical to that shown in fig. 1. For massless continuum fermions, 
the integrals ( [A 3D may be evaluated analytically and yield 



C 



18 



T 



D(t) 



3 



sinh(27rTt) 



(2vrTt + ztt) (coth 2 (27rTt) - §) - coth(27rTt)j . (A4) 



APPENDIX B: LONG-TIME DECAY OF D(T) 

The long-time behavior of D{t) may be extracted by applying a stationary-phase approx- 
imation to the three-dimensional momentum integral ( |2.4cj ) defining D(t), or equivalently by 
analyzing the effect on the one- dimensional Fourier transform (|2.7|) of the non-analyticities 
in the spectral density p vv (uj) ( |2.<j| ). In the massive continuum theory, the only saddle point 
of e 2luJpt occurs at p = and the resulting saddle point contribution is 

D(t\ ~ 2 o„ / 2i[m+p 2 /2m+0(p 4 /m 3 )]t 

a J (27r) 3 m z m 2 

T 2 t~ 5 / 2 / / 1 \\ 

cos(2mt-|7r) (l + O (— . (Bl) 



(47rm) 3 / 2 4 V Vwit 

At intermediate times, (3 <^it m~ 1 , one may show that 

T 2 m 

D(t) cont = D(t) massless 

(1 + 0{m 2 (3 2 )) - — (1 + 0(m 2 t 2 )) , (B2) 

where -D(t) mass i ess ~ -^T 3 e~ 47TTt (2rrtT — 1 + m) is the massless continuum result of (|2.9|) . 

For the lattice theory (on a simple cubic lattice), the dispersion relation uj p has additional 
saddle points at p = (0, 0, ±7r), (0, ±7r, ±7r), and (±7r, ±7r, ±7r), on the edges and corners of 
the Brillouin zone, with frequencies of 2/a, 2\/2/a, and 2y / 3/a, respectively. Each saddle 
point produces an oscillating t~ 5 / 2 contribution at the associated frequency, so that 



T 2 a 3 / 2 

D(t) lattice ~ D(t) cont H ^ 5/2 



cos(4t/a + f 7r) cos(4\/2 t/a - |7r) cos(4\/3t/a + |7r) 
(8vr) 3 / 2 + (8V2n) 3 / 2 + (Sv^vr) 3 / 2 
x (1 + O (a/t)) , (B3) 



where -D(t) cont behaves as shown in ( [B 1| ) for £ 3> m _1 , or the second term of ( |B"2| ) for 
a In <C £ <C mT x . Note that, for times not much larger than a, the amplitudes of the lattice 
oscillations in ( |B"3| ) are comparable to the 0(T 2 /a) steady-state value C of the (perturba- 
tive result) for the current correlator. At times of order 1/m, where the continuum mass 
oscillations begin to be evident, the amplitude of the oscillation has dropped to 0(T 2 m), or 
0(ma) = 0(/3 L ) times the steady-state value. 
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